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Charged soft-matter systems-such as coUoidal dispersions and charged polymers-are dominated 
by attractive forces between constituent like-charged particles when neutralizing counterions of 
high charge valency are introduced. Such counter-intuitive effects indicate strong electrostatic 
coupling between like-charged particles, which essentially results from electrostatic correlations 
among counterions residing near particle surfaces. In this paper, the attraction mechanism and 
the structure of counterionic correlations are discussed in the limit of strong coupling based on 
recent numerical and analytical investigations and for various geometries (planar, spherical and 
cylindrical) of charged objects. 
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I. INTRODUCTION 



Electrostatic interactions and processes involving electric charges appear ubiquitously in biological and soft-matter 
systems. Electric charges make materials water-soluble and lead to many important technological and biological 
applications. In industry, for example, the stability of colloidal dispersions is often a desirable property as in the 
case of paints and food emulsions such as milk. One way to stabilize colloidal suspensions against coagulation or 
flocculation (occurring due to attractive van-der-Waals forces) is to generate long-range repulsive interactions between 
constituent colloidal particles by imparting permanent like-charges to these particles Q, |^ |^ , or by grafting charged 
polymeric chains to their surfaces (forming hairy particles or polymeric brushes) 0. Charged polymers, or the so- 
called polyelectrolytes, and their synthesis have also attracted a lot of attention [3 since, due to their water-solubility, 
they offer a useful option in design and processing of non-toxic environmentally-friendly materials. In biology, on 
the other hand, electrostatic effects emerge in many striking examples such as the DNA-packaging process in the 
cell nucleus. In each human cell, a total length of about 2m of DNA-which bears a total negative charge of about 
10^"e {i.e. one elementary charge, e, per 1.7A)-is stored inside the cell nucleus with a diameter of less than lO/im. 
This storage process involves a hierarchical structure on the lowest level of which, short segments of DNA are tightly 
wrapped around positively charged proteins of few- nanometer size (the so-called Histones). It is shown experimentally 
[Pi that such a tightly wrapped state is only stable for intermediate, physiological salt concentrations, at which an 
optimal balance between self-repulsion of DNA segments and the DNA-Histone attraction is achieved. 

Charged macromolecules (macroions) in solution are always surrounded by neutralizing counterions, and also in 
general by colons. Counterions form electrostatically-bound clouds in the proximity of macroions and in many cases, 
predominantly determine the electric properties of charged solutions In particular, counterions can alter the 



effective interaction between like-charged macroions 
them in certain physical conditions [3,[ 



■enerate a dominant electrostatic attraction between 
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Like-charge attraction manifests itself in a number of famous 
examples, namely, the condensation of DNA molecules bundle formation of stiff polymers and aggregation 
of colloidal particles |3^ [s^, 13- Interestingly, such an attraction emerges only in strongly-coupled systems, e.g. 
when macroions are highly charged (with surface charge densities up to le/nm^ as in the DNA system), neutralizing 
counterions are multivalent, or when the temperature or the dielectric constant of medium is low. For instance, the 



2 



DNA condensation process, in which long DNA molecules condense into a tightly packed, circumferentially wound 
torus, is observed in experiments where multivalent counterions (such as trivalent spermidine ions) are introduced 
[T?il |. A similar trend has also been found in numerous numerical simulations of like-charged membranes, colloids and 
polymers ©HIlllilllllMlillillillilliillSlllllilllH where highly 

charged macroions are found to form closely-packed bound states due to attractive forces of electrostatic origin. 
These attractive forces are of typically large strength compared to the usual van-der-Waals attraction and may have 
significant practical implications where, for instance, multivalent counterions are present. 

In the weak-coupling condition, e.g. when surface charge densities and charge valency of counterions are small, 
macroions are found to repel each other. In fact, from a theoretical point of view, one can argue that weakly-coupled 
systems should be described well by the mean-field approximation. Therefore, the effective interaction between 
macroions may be calculated using the mean-field solution for the electrostatic potential field in space, which follows 
from the so-called Poisson-Boltzmann (PB) equation that has been widely studied in the past It is known that 

the Poisson-Boltzmann theory-despite its success in many applications-can only yield repulsion between like-charged 
macroions [tUsoIIsHI. 

The main scenarios which are put forward to explain the phenomenon of like-charge attraction have gone beyond the 
mean- field fovel^by^m^n^ j^^ ^*^^ 

correlations are systematically neglected on the mean-field level but become enhanced in strong-coupling conditions 
mentioned above. Recent theoretical attempts to incorporate ionic correlations include integral-equation methods 
I25l l4p^, p erturbative improvement of the mean-field theory including Gaussian-fiuctuations theories 0, lisl l46l l47l 
l48l l49l Isol Isil Is^ Is^ . IsS. IsM Is^ Is^ . and local density functional theory 59, 60], which compare well with numerical 
simulations and all exhibit attraction. These methods are mostly applicable for large separations between macroions 
or in the regime of low coupling strengths (the so-called high-temperature regime), and can not characterize the 
closely-packed bound state between like-charged macroions. An alternative approach was triggered by Rouzina and 
Bloomfield 6l] with the insight that counterions form two-dimensional highly-correlated layers at macroionic surfaces 
for high coupling strengths (the so-called low-temperature regime). Such structural correlations give rise to attractive 
forces [62, EE EE EE EE EE IZE 1^ j which mainly result from ener getic origins and can account for the 

closely-packed bound state of like-charged macroions HE IM EE EE EE El EEH EE EE EE EE Iz3 • 

In this paper, we shall chiefly consider this latter regime of strong coupling. First we shall briefly discuss the 
regime of parameters where electrostatic correlations are expected to be small or large in terms of the main length 
scales in a charged system (Section In Sections IIIIIVI we shall focus on the system of counterions at one or two 
planar charged walls to demonstrate the gross physical picture in both weak-coupling (or mean-field) and strong- 
coupling regimes. The role of curvature and its main physical consequences are considered in Sections IVIIVIIll As 
we shall see, the so-called counterion-condensation process occurring at curved surfaces plays an important role in 
determining the attraction regime between spherical (Section rVIII) and cylindrical (Section^TOJ macroions. We focus 
on analytical results only based on two exact asymptotic theories, namely, the mean-field Poisson-Boltzmann theory 
and the strong-coupling (SC) theory. The crossover regime between these two limits are considered using recent 
numerical simulations, which are useful to determine the regime of applicability of both asymptotic theories, and also 
to examine the validity of their predictions. 

We restrict our discussion to a primitive model in which particles only interact with Coulombic forces and in some 
cases, also with short-range excluded-volume repulsions. The inhomogeneous charge distribution of macroions and 
also ion-surface adsorption effects [t^I are neglected. We do not consider colons and take the solvent effects into 
account only through the continuum dielectric constant of medium e (for water at room temperature e « 80), which is 
mainly assumed to be uniform and equal everywhere in space. The role of the dielectric jump at charged boundaries 
will be discussed briefly in Section for the case of two charged walls. 

II. LENGTH SCALES IN CHARGED SYSTEMS: FROM MEAN-FIELD TO STRONG-COUPLING 

REGIME 

To distinguish the regimes of parameters where attractive or repulsive forces may arise between macroions, one 
needs first to study the length scales that appear in a classical charged system. Let us consider a system of macroions 
with uniform surface charge density of —as (in units of the elementary charge e) and neutralizing counterions of 
charge valency q at temperature T . (Hereafter, we conventionally assume that macroions are negatively charged and 
counterions are positively charged, thus Us and q are both positive by definition.) 

A characteristic length scale in such a system is set by comparing the thermal energy scale, k-^T^ with the Coulombic 
interaction energy between counterions, V{r) — q^e^ / (47reeor), where r is the distance between two given counterions. 
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FIG. 1: a) For small coupling parameter, S ^ 1, counterions form a diffuse three-dimensional layer, b) For large coupling 
parameter, H ^ 1, the counterionic layer is essentially two dimensional since the typical lateral distance between counterions, 
a±, becomes much larger than the Gouy-Chapman length, fi. In this regime, counterions are strongly correlated and surrounded 
by a correlation hole of size ~ ax ^. 



The ratio between these two quantities may be written as V/{kBT) = q^i-g/r, where 

is the so-called Bjerrum length, which measures the distance at which two elementary charges interact with thermal 
energy k^T (in water and at room temperature €b ~ 7.lA). Thus the rescaled Bjerrum length 

4 = <?'^B (2) 

may be taken as the relevant length scale to characterize the strength of mutual counterionic repulsions against 
thermal fluctuations in the system. 

Other length scales are set by considering the charge distribution and the specific geometry of macroions. For 
simplicity, let us concentrate here on a planar system composed of a planar charged wall of infinitely-large extension 
and neutralizing counterions confined to one half-space (Figure^. This model is physically relevant for charged 
membranes and also for macroions with large radii of curvature, which behave like plates at small distances from their 
surface. 

In this system, an additional length scale may be obtained by comparing the thermal energy k^T with the energy 
scale of the counterion-wall attraction, U{z) = (7(Tse^z/(2e£o)i where z is the vertical distance from the wall. Hence 
we have the ratio U/{kBT) — z/fj,, where 

gives the so-called Gouy-Chapman length, which measures the distance at which the thermal energy equals the 
counterion-wall interaction energy. The Gouy-Chapman length also gives a measure of the thickness of the counterionic 
layer at a charged wall, (z) ~ fi, as we shall see later. 

In principle, one may tune the system parameters in such a way that either of the two length scales £b and fi 
become arbitrarily large or small. In fact, only the dimensionless ratio between these two quantities is relevant and 
uniquely describes the physical regimes of the system l82l| , namely, 

S = - = 2V4^s, (4) 

which is known as the electrostatic coupling parameter. For small coupling parameter S ^ 1, the counterion-wall 
system has a relatively large Gouy-Chapman length (or small Bjerrum length), which indicates a loosely bound 
counterion cloud at the charged wall (Figure^). For large coupling parameter S ^ 1, the Gouy-Chapman length is 
relatively small (or the Bjerrum length is large) indicating that counterions are strongly attracted toward the charged 
wall (Figure ^d). 

Further insight into the structure of the counterionic layer may be obtained by considering the typical distance 
between counterions at the surface. For counterions residing near the charged surface, the local electroneutrality 
condition implies a typical lateral separation of 

a± ^ W — , (5) 

V O-s 
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TABLE I: Typical values of physical parameters for realistic charged systems: as and R denote the surface charge density 
and the radius of curvature of charged objects, q is the charge valency of counterions, and fi, H and ^ are the Gouy-Chapman 
length, fi = 1/ {2-Kql-B<Js) (Eq. @), the coupling parameter, H — q^^b/m (Eq. ||1J), and the Manning parameter, — R/ ^ 
(Eq. 1261 '). respectively. (The role of curvature and Manning parameter for cylindrical and spherical macroions is discussed in 
Sections I VII Villi ) The Bjerrum length is taken here as 1-b ~ 7.lA corresponding to an aqueous medium of dielectric constant 
e = 80 at room temperature. 



since each counterion neutralizes the charge of an area given by a\ ^ q/ds (up to a geometrical prefactor of the order 
one). Comparing this length scale with the Gouy-Chapman length, we have 



(6) 



Hence in the strong-coupling regime, S 3> 1, counterions essentially form a quasi two-dimensional layer as their lateral 
separation at surface becomes much larger than the Gouy-Chapman length a± ^ fi (Figure^). On the other hand, 
the structure of such a layer is dominated by mutual repulsions between counterions, which freeze out lateral degrees 
of freedom. Hence counterions become laterally correlated and surrounded by a large correlation hole of size a± from 
which neighboring counterions are statistically depleted [6ll 163. l66l| (see Section IHI CI for a more detailed analysis 
|39j). This indicates a trend toward crystallization in the ionic structure for increasing coupling parameter, which can 
be corroborated by considering the effective plasma parameter relevant for this situation 



r = — =1/2 



(7) 



The parameter F gives a measure of mutual Coulombic repulsions between counterions at a neutralizing surface (the 
2D one-component plasma) [6ll l64l l6^ . For increasing S, F increases while the average counterion- wall interaction 
(per fceT) remains of the order of unity, {U)/{k'QT) ^ {z) I ijl ~ 1. The Wigner crystallization of the 2D one- 
component plasma is known to occur for F > Fc ~ 125 |83| . which corresponds to the range of coupling parameters 
S > « 3.1 X 10^ m (see Section UlTnil . 

In the weak-coupling regime (S <^ 1), no crystallization is expected to occur and the counterion layer has a 
three-dimensional fluid- like structure (ax ^ n) |84| (Figure^). Thus the two asymptotic regimes of weak coupling 
(S ^ 1) and strong coupling (S ^ 1) may be distinguished physically by the structure of counterionic layers at 
charged surfaces. In Sections IIIII and HVl we shall briefly review the main results obtained in each of these regimes for 
the classical example of counterions at one and two charged walls (planar double layers) . 

But before proceeding further, it is useful to consider the typical values of the coupling parameter in realistic 
systems. In Tabled we show few typical examples of both weakly-coupled and strongly-coupled systems. As already 
seen from Eq. the coupling strength grows quite rapidly with the counterion valency (S g^), which agrees with 
experimental and numerical evidence indicating highly growing correlation effects for increasing counterion valency 

SIpLRMEIillillillillillilE 
liol l4lL |42 | . In fact as known from these studies, typical coupling strength of 2 ~ 10^ (or larger) already reflects 
strong-coupling regime and a value of S ~ 1 (or smaller) typically corresponds to the weak-coupling regime. 
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III. COUNTERIONS AT A CHARGED WALL 
A. Weak-coupling (or Mean-field) regime H <C 1 



For small coupling strength, one may employ the mean-field approximation to describe the counterionic layer since 
each comiterion interacts with a diffuse cloud of other counterions. The mean-field ap pro ximation systematically 
neglects correlations among counterions and is formally valid in the strict limit of S — > |55|. 

The mean-field theory is governed by the so-called Poisson-Boltzmann (PB) equation QHSi Issf 



££0 



exp{-qeip/kBT), 



(8) 



which is to be solved for the mean electrostatic potential in space, ip, using proper boundary conditions at macroion 
surfaces. The corresponding density profile of counterions, /3pb(x), follows from the solution of the PB equation (jS)) 
and using the relation /9pb(x) = poexjp{—qe'ip/kBT), where po is a normalization prefactor. 

For the system of point-like counterions at a single uniformly-charged wall (in the absence of salt) , the PB theory 
predicts an algebraically-decaying density profile of the form 0, Ql 



1 



Ppb{z) ^ 

27r^BCT2 {z/n + 1)2 



(9) 



where z is the distance from the wall. The density of counterions at contact is obtained as p{z = 0) = 27r£BO's , which 
is an exact result within the present model and valid beyond the mean-field level . As seen the PB theory predicts 
an extended counterionic density profile (with diverging moments) in agreement with the qualitative considerations in 
Section|n]for a weakly-coupled system. Nonetheless, the PB density profile ^ is normalizable to the total number of 
counterions reflecting the fact a charged wall binds all its counterions. Note also that the Gouy-Chapman length, p,, 
in this case equals the height of a layer at the wall which contains half of the counterions, and thus may be associated 
with the typical layer thickness at low couplings. 



B. Strong-coupling regime ^ 2> 1 

In this regime, the liquid-like ordering (or crystallization at sufficiently large couplings) of counterions renders the 
mean-field theory an invalid des crip tion of the system. Yet one can obtain a simple analytical description for the 
counterionic layer as follows |^ jg^ • 

Since for S ^ 1 counterions become highly separated from each other in a quasi-2D layer at the charged wall 
(Figure^), one may consider the system as a collection of laterally frozen correlation cells, each consisting of a single 
counterion interacting with an area of the wall of size ~ a± (Eq. (O)- Since in this regime a±_ 3> p, the dominant 
contribution to the density profile of counterions at the wall is obtained by considering only the vertical degree of 
freedom, z, through which single counterions are coupled to the wall with the interaction potential U/{k-BT) w z/p. 
Hence using the Boltzmann weight, one has the following density profile 

Psc{z) = Poexp{-z/p). (10) 

The prefactor in the above expression (the contact density) may be fixed from the normalization condition for the 
density profile and the global electroneutrality of the system as po = 27r^BO's . The strong-coupling density, psciz), 
drops quite rapidly as one moves away from the charged wall. The average distance of counterions is obtained to be 
equal to the Gouy-Chapman length, (z)sc = A*- 

The above density profile, which essentially follows from single-particle contributions, was obtained by Shklovskii 
[sS. l6^ using a Wigner-crystal model for large Coulombic coupling. The asymptotic analysis of Ref. showed 
that the partition function of the system for S — > cxd adopts a series expansion in powers of 1/S, the leading term of 
which is given only by single-particle contributions. The multi-particle contributions enter in higher-order terms (in 
the form of a virial expansion). The leading term defines the asymptotic strong-coupling (SC) theory, which for the 
counterion-wall system gives exactly the density profile H10|l . 



C. Intermediate-coupling regime 



In realistic systems, the coupling parameter is always finite (see Table QJ. Therefore, it is important to examine 
whether and how the preceding analytical results for the two limits of strong (S — > oo) and weak (S — > 0) coupling 
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FIG. 2: Density profile of counterions at a cfiarged wall as a function of the distance from the wall, z. The symbols are the 
data from Monte-Carlo simulations for H = 0.1 (open diamonds), H = 10 (open triangles) and H = 10* (open stars). The solid 
curve and the dashed curve represent the predictions of the mean- field (PB) theory, Eq. @, and the strong-coupling theory, 
Eq. lion , respectively. The density profile is rescaled by its values at contact p(z = 0) = 27r£BO"s a^nd the distance from the 
wall is shown in units of the Gouy-Chapman length, [i. 



may be applied to such systems. A useful approach to investigate the regime of intermediate couplings is to employ 
numerical simulation methods. Other methods include systematic improvement of both mean-field and strong-coupling 
theories that will be considered only for the interaction between two walls in Section IWI 

The density profile of counterions at an infinite ly l arge charged wall has been calculated for various coupling 
parameters using Monte-Carlo simulations in Ref. |37|. (The simulation model is similar to what we described in 
Section ITTl (Figure ^1, where counterions are taken as p oint-like particles confined to one half-space and periodic 
boundary conditions are used in lateral directions [39j|.) In Figure [3 we show the simulated density profile of 
counterions (symbols) for coupling parameters S = 0.1,10 and 10* along with the analytical predictions of the 
mean-field PB theory (sohd curve, Eq. Q) and the strong-couphng theory (dashed curve, Eq. (fTUIl ). Note that in the 
Figure the distance from the wall is rescaled with the Gouy-Chapman length; hence, the PB and SC profiles appear 
to have a similar decay length of about unity in rescaled units. As seen, the data crossover from the mean-field (PB) 
prediction to the SC result by increasing S, quantitatively confirming the validity of both theories at small couplings 
(about S — 0.1) and large couplings (about S — 10*) respectively. 

Further analysis of the simulated density profile |33 reveals a distance-dependent crossover at intermediate coupling 
parameters: while the large distance behavior of the system is described well by the mean-field theory, at sufficiently 
small distances from the wall, the system roughly follows the strong-coupling prediction at intermediate S. This is 
in fact a quite general property that we shall investigate in more detail for the interaction between two charged walls 
in Section HVl The crossover in the single- wall system can be understood qualitatively by noting that the asymptotic 
strong-cou pling density profile, Eq. IjlOl) . remains valid only within distances smaller than the correlation hole size 
z < aj_ |39ll64L l69l l86l| . Thus at finite coupling parameters, the strong-coupling regime may be characterized by 

i<sV2, (11) 

where we have used Eq. (O . For distances, z, from the wall which are larger than a_L, lateral interactions between 
counterions become relevant and modify the density profile. For z ^ , one intuitively expect that dominant many- 
body effects lead to a mean-field-like behavior as shown in previous studies 1391 |69 l . In fact, the mean-field results 
turn out to be valid for distances > S (up to some logarithmic corrections) 'SJI, '6!l|. 

An interesting problem is to examine the development of correlations between counterions (including the formation 
and size of the correlation hole) as the coupling parameter increases. This goal may be achieved by considering the 
two-dimensional pair distribution function of counterions, which is defined through 

52D = ^l^^K^^xy-r^^y^. + v^y^^y), (12) 

where the sum runs over pairs of particles, v-j^y is a 2D vector and 'ixy,i is the lateral position of particle i in the 
xy-plane (see Figure^. Physically 1720 gives the ratio between the probability of finding two counterions at distance 
Txy = \vxy\ from cach other and the expected probability for a homogeneous 2D fiuid with the same bulk density. 
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FIG. 3: a) The two-dimensional pair distribution function of counterions at a charged planar wall plotted as a function of the 
lateral distance between counterions (the distribution function is obtained by averaging over the height z). The symbols are 
the data from Monte-Carlo simulations for H = 1 (filled triangles), H = 10 (open triangles), H = 100 (filled diamonds) and 
H = 10* (open stars). The lateral distance is rescaled by the length scale 2a^ — 2^ (//(ttcts), which gives a measure of lateral 
separation between counterions (see the text). 

b) The simulated excess heat capacity of the system of counterions at a charged wall (per number of counterions, A^) plotted as 
a function of the coupling parameter, H. (The number of counterions in this case is A*' = 100 in a square simulation box, which 
is periodically replicated in lateral directions.) The broad hump at intermediate couplings (10 < H < 100) reflects the structural 
change in the counterionic layer due to increasing correlations between counterions. At large coupling strength (H « 3.1 x 10'*), 
the counterionic layer undergoes a crystallization process indicated by a pronounced peak in the heat capacity. 



The Monte-Carlo results for this quantity are shown in Figure |39j . It shows that for small coupling parameter 
(S = 1, filled triangles), there is only a very short-range depletion zone at small separations between counterions. But 
a pronounced correlation hole is created at the regime of coupling parameters 10 < S < 100, where the distribution 
function vanishes over a finite range at small inter-particle separations. For larger coupling strengths, the correlation 
hole becomes more pronounced and is followed by an oscillatory behavior in the pair distribution function (S = 10^, 
open stars) . This indicates a liquid- like order in the counterionic structure in agreement with qualitative considerations 
in Section m Note that the distance coordinate in Figure |2t is rescaled with 2a° , where a*[ = ^q/{iT(Js) is obtained 
by assuming that the area of the wall neutralized by a counterion has a circular form (Figure^). The location of the 
first peak of (72D for S = 10"* appears at a distance of rxy/{2a!\^ k, 0.9, which roughly gives the typical lateral distance 

between counterions. In a perfect hexagonal crystal, the peak is expected to occur at r^y/ (2a*]^) = \J 7r/(2\/3) ~ 0.95, 

and in a perfect square crystal at rxy/{2a!\^ = VtF/^ ~ 0.89. The crystallization is in fact found at even larger 
coupling parameters |87j |. which may be seen also from the behavior of the heat capacity for increasing S. 

In Figure l^b, the simulated excess heat capacity of the counterion-wall system (obtained by omitting the trivial 
kinetic energy contribution 3A;b/2) is shown for various coupling parameters. The crystallization of counterions at 
the wall is refiected by a pronounced peak at large coupling parameters about « 31000, which roughly agrees 
with the Winger-crystallization threshold of a 2D one-component plasma [s^- The characteristic properties of the 
crystallization transition in the counterion-wall system are yet to be specified, which requires a detailed finite-size 
scaling analysis in the vicinity of the transition point. Another interesting behavior is observed in Figure at the 
range of coupling parameters 10 < S < 100, where the heat capacity exhibits a broad hump. This hump does not 
represent a phase transition (3^, but it is associated with the onset of the correlation hole around counterions and 
the structural change in the counterionic layer from being a three-dimensional layer at low couplings to a quasi-2D 
layer at large couplings. In the region between the hump and the crystallization peak (for 200 < S < 10''), the heat 
capacity is found to increase almost logarithmically with S. The reason for this behavior is at present not clear. 



IV. REPULSIVE AND ATTRACTIVE LIKE-CHARGE INTERACTIONS: PLANAR GEOMETRY 

Now let us consider the interaction between macroions in the two regimes of mean field and strong coupling. In 
this Section, we focus on the the planar system of two parallel charged walls of uniform surface charge density — (Tg 
at separation A from each, where q-valent counterions fill only the space between the walls (the dielectric constant 
is also assumed to be uniform in space)-see Figure 0] This model is relevant for interaction between two charged 
membranes or two macroions of large radii of curvature. 

In this system, an extra length scale is set by the wall separation, A. Two limiting regimes of repulsion and 
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FIG. 4: Schematic representation of the asymptotic interaction regimes of a) mean field and b) strong couphng for two hke- 
charged waUs. The mean-field regime is obtained at large separations between the walls (compared to other length scales) and 
is dominated by the repulsive osmotic pressure of counterions. For small wall separation (compared with the typical counterion 
spacing), the walls attract each other since counterions are isolated in correlation cells of large lateral extension ~ a_L/2 2> A 
(shown by a dotted loop) and mediate a dominant single-particle attraction between the walls. 



attraction may be distinguished qualitatively by comparing A with other length scales of the system as follows. 



A. Mean-field regime: repulsion 



First consider the limit where the wall separation, A, is large compared with all other length scales in the system 
and also that the system is weakly coupled, S ^ 1 (see Figure 0^). In this case, counterions form a diffuse layer 
at each wall, but due to large wall separation, the system is approximately decoupled into two nearly neutral sub- 
systems, each consisting of a charged wall and its counterionic cloud. The effective pressure acting between the walls 
is dominated by the osmotic pressure of counterions across the mid-plane, since the overall electrostatic field at the 
mid-plane is zero due to the charge neutrality of each sub-system. This osmotic pressure is positive and therefore 
gives an effective repulsion between the walls. 

The mid-plane osmotic pressure is proportional to the local density of counterions, /Omid, following the ideal-gas 
equation P ~ PmidkBT, where pmid drops roughly with the inverse square of the wall separation for large A as it 
follows from Eq. I^. This yields the scaling form of the repulsive pressure between the walls as P(A) ^ A~^. The 
formal derivation of the pressure based on the PB equation supports the above result for large separation. The full 
PB solution for arbitrary A follows as 0,0,113 



/3Ppb(A) 



A (13) 



(with (3 = I/ZshT), where A is determined from the transcendental equation A^/^ tan[A^/^(A/2/x)] = 1. For large 
A/fi^l, the PB solution yields 



(14) 

which is expectedly independent of the surface charge density of the walls. 



B. Strong-coupling regime: attraction 



Now we consider a different asymptotic regime in which the system is strongly coupled, S ^ 1, and the distance 
between the walls is smaller than the lateral spacing between counterions at each wall a± ^ A (Figure^). Since 
counterions are highly separated from each other, the two opposite layers of counterions tend to form an inter-locking 
pattern at equilibrium, where each counterion from one layer faces a bare area of the other wall with nearly equal 
(but opposite) charge. It is evident that this pattern leads to an effective attractive force between the walls with a 
purely energetic origin. The effective pressure acting between the walls may be estimated using a simple argument 
that takes into account the highly correlated structure of counterions 

In the asymptotic limit considered here (i.e. for S 3> 1 and a± ^ A), each counterion is essentially confined 
and isolated in a "correlation cell" that consists of a single counterion sandwiched between two opposing sections 
of the walls with lateral size of about a±/2 (Figure 21d). Since a± ^ A, the effective pressure between the walls is 
dominated by the contribution coming from each single correlation cell and lateral interactions between these cells 
may be neglected. 
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FIG. 5: a) Effective pressure between two like-charged walls as a function of their distance, A. The symbols are MC simulation 
data for H = 0.5 (open diamonds), H = 10 (filled diamonds), H = 20 (open squares), H = 100 (open stars) and H — 10^ (open 
triangles). The solid curve is the prediction of the mean- field PB theory, Eq. 113II . and the dashed curve is the strong-coupling 
prediction, Eq. 1)15^ . The pressure and the wall separation are shown in rescaled units as indicated on the graph, 
b) Attraction and repulsion regimes shown in terms of the rescaled wall distance, A/fi, and the coupling parameter, H. Symbols 
are the simulation results indicating the zero-pressure points (connecting lines are guides to the eye). Filled symbols show the 
thermodynamically stable bound state of the two walls at small separations. Open symbols indicate the meta-stable or unstable 
states of zero pressure. Attraction sets in for H > 12 and a first-order phase transition occurs at H « 17. 



The electrostatic energy of the system per cell is the sum of the interactions between the two surfaces with each 
other and with the single counterion, which-using the electroneutrality condition and the fact that the wall separation 
is small-follows as /3ue ~ 27r£BO's A per unit area. This energetic contribution gives an attractive pressure of /3Pe ~ 
— 27r^BO's between the walls. On the other hand, the entropic contribution due to the counterion confinement is of 
the order S ^ fee In A (per cell), which generates a repulsive component. The total pressure between the walls is then 
obtained by combining these two effects and may be written as 

/3^sc(A) 2^, 

This expression clearly predicts a closely-packed bound state between the walls with equilibrium surface separation, 
A*, being equal to twice the Gouy-Chapman length, i.e. 

A* 2^. (16) 

The like-charged walls attract each other for A > A, and repel at smaller distances. 

The energetic attraction in the strong-coupling regime, which results from structural correlations, was first obtained 
by Rouzina and Bloomfield and investigated later by several workers lilililiilillelillillTiEi- The 
rigorous derivation of the expression H15|) for the pressure was given using the asymptotic strong-coupling theory |69j| , 
which shows that Eq. H15|l is the exact result in the limit of S ^ cx). 

Note that in the presence of counterions of finite diameter ad, the equilibrium separation between the walls, Eq. 
(|16|l . increases by an amount equal to the counterion diameter and reads 

A, = CTci + 2^ (17) 

assuming that counterions interact with a hard-core excluded- volume interaction with the walls. Clearly, the excluded- 
volume interaction between counterions themselves is irrelevant in the strong-coupling regime since counterions are 
highly separated from each other {a± ^ A). 



C. Numerical simulations and the crossover regime 

In order to examine the preceding asymptotic results in the mean-field and strong-coupling regimes, we consider 
the Monte-Carlo simulations of the two-wall system, which enable one to investigate the mechanism of like-charge 
interaction beyond the above limiting cases. 

FigureE^ shows the simulated effective pressure acting between two like-charged walls (in the presence of point-like 
counterions) for various coupling parameters |38l l39j. The pressure becomes negative, and thus indicates attraction 
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between the walls at small and intermediate separations, when the coupling parameter exceeds an intermediate 
threshold (see below). The onset of attraction at intermediate couplings (10 < S < 100) agrees with the onset of 
correlations between counterions as discussed in Section lTlI CI For small couplings (S = 0.5, open diamonds), the data 
quantitatively support the mean- field PB prediction (solid curve), Eq. (|13|l . and for very large couplings (S = 10^, 
open triangles), they agree with the SC prediction (dashed curve), Eq. (|15|l . quite well. 

The behavior of the pressure may be summarized in a phase diagram as shown in Figure Eh, which shows the 
regions of positive (repulsive) and negative (attractive) pressure separated by a line of zero pressure. The attractive 
region only appears for coupling parameters S > 12. The filled symbols (connected with a solid line) show the stable 
bound state of the two walls, while the open symbols (connected with a dashed line) correspond to meta-stable or 
unstable states of the two wall {i.e. the local minima or the maximum of the free energy of the system, where the 
free energy is obtained from the data by integrating the pressure from infinite distance to a finite distance, A) |39l |. 
As seen, the stable bound state exhibits an equilibrium wall separation quite close to the strong-coupling prediction, 
A*//i = 2 (Eq. (|16|l ), for moderate to large coupling parameters. The thermodynamic behavior of this system has 
been studied in Ref. 1 which predicts a first-order unbinding transition at S « 17. 

The effective pressure also exhibits a distance-dependent crossover at intermediate couplings ^3§, 's^, '6^ : at small 
wall separations, the data closely follow the SC curve, while for large separations, they tend to the PB curve and 
display a mean-field-like repulsion. In order to study the crossover behavior analytically, one needs to consider the 
extension of both asymptotic theories of mean field (S — )■ 0) and strong coupling (S — > oo) to finite-coupling situations. 



1. loop expansion: sub-leading corrections to the mean-field PB theory 

The mean-field PB theory is obtained from a saddle-point approximation in the limit of S ^ [s^. Therefore, 
one way to incorporate finite-coupling effects on a systematic level is to calculate the higher-order corrections to the 
saddle-point solution by means of a loop expansion. The loop parameter turns out to be the coupling parameter S 
and the effective pressure between the walls may be expanded about the mean-field solution as 

P(A) = Ppb(A) + EPj^'^iA) + O(S^), (18) 
where PpB is the PB solution (|13|l , and Ppg is the first-loop or the Gaussian correction term Q3 



27r£Bcr? VA 



C(3) 



4 4 



— +7:' ln(A/^^) 



(19) 



Clearly, the Gaussian correction term contributes an attractive component, which comes from correlations between 
fluctuations in the counterionic clouds at opposite walls. These fluctuations tend to polarize each other giving rise 
to attraction in the same way as other fluctuation-induced attractive forces (such as dispersion interactions) are 
generated [54| . It is tempting to argue that the Gaussian correction term turns the net pressure between the walls 
into an attractive pressure for large enough S. However, the onset of attraction in fact signals the break-down of 
the loop-expansion scheme as used above, since the correction term becomes comparable to the leading PB term 
|69|| . Therefore, the Gaussian-fluctuations picture remains valid only at sufficiently small couplings (or the so-called 
high-temperature regime) and also for sufficiently large separations A//i 3> 1 (see the discussion in Refs. |5^l68l 69]). 

The regime of validity of the loop expansion (and that of the mean-field PB theory) at a finite coupling parameter, 
S, may be estimated by comparing the sub- leading and the leading terms in Eq. (|18|l . that gives ,,39> i69j 

t4P^ > (20) 
ln(A/^) ^ ' 



2. virial expansion: sub-leading corrections to the strong- coupling theory 



In the strong-coupling regime, the finite-coupling corrections may be taken into account using a virial-expansion 
scheme, which is obtained as a series expansion in powers of 1/S about the asymptotic strong-coupling solution (for 
S ^ oo) |69| . The effective pressure between two like-charged walls adopts the following large coupling expansion 

P(A) = Psc(A) + ip«(A) +0(S-2), (21) 



11 




FIG. 6: a) Regimes of applicability of the asymptotic theories of mean field Poisson-Boltzmann and strong coupling for the 
system of two like-charged walls and counterions. There is an intermediate regime of rescaled distances between the walls, 
A//i, and coupling parameters, H, where finite-coupling effects can not be captured by series expansions around the mean-field 
or the strong-coupling solutions. Possible physical situations in this regime are schematically shown in b) and c)-see the text. 



where Psc is the SC prediction, Eq. ifT^ . and P^q is the first correction term [sMIs^IsS 



(1) 

■ = — (22) 

which contributes a repulsive component to the total pressure. This finitc-couphng correction can be used only at 
sufhciently large couplings and small wall separations where the correction term itself is small, i.e. where the 1/S- 
expansion scheme remains valid. One may estimate the regime of validity of this expansion (and thus the regime of 
applicability of the SC theory) at a finite coupling parameter from Eqs. and (|^ as ,£21 



< E. (23) 



This estimate in fact agrees with our qualitative discussion in Section IIVBI which predicts the strong-coupling 
attraction for small wall separation compared with the lateral distance of counterions, a±, that is for 

A < a^, (24) 

which in units of the Gouy-Chapman length (and using Eq. ^) reproduces Eq. H23|l . Note also that the equilibrium 
wall separation predicted by the SC theory, A*//^ — 2 (Eq. fulfills the above criterion for S > 4. 

Equation H24ll -or H23|l in rescaled units-is also known as the Rouzina- Bloom field criterion 61], which is established as 
a generic attraction criterion for highly-charged macroions including charged spheres and cylinders [30ll33. l33L l4ll FH . 

The above discussions may be summarized in a diagram as shown in Figure IHt^ specifying the range of parameters 
(coupling parameter and the wall separation) where the strong-coupling or the mean-field picture prevails. As seen 
there appears a gap in the diagram, where neither of the theories can be extended to include finite-coupling effects 
via the series-expansion methods mentioned before. The physical situations to which this gap corresponds have been 
illustrated in Figures IHh and E); (compare these Figures with Figure 0}. Figure (Sja shows a system in which the 
Gouy-Chapman length is the smallest length scale and the wall separation is large such that /i < a_L < A. In units 
of the Gouy-Chapman length, we have 1 < S < (A//i)^. In this case, the PB approach is not valid and in a rough 
approximation, the two layers are decoupled and each layer is separately described by the strong-coupling density 
profile for a single wall. Yet a systematic theory for the effective interaction in this regime is missing. Figure EJ; 
shows a system in which the Gouy-Chapman length is the largest length scale and < A < /i, or in units of the 
Gouy-Chapman length, S < (A//i)^ < 1. In this case, counterions form a confined gas with local three-dimensional 
correlations for finite S. Interestingly in this regime both SC theory and PB theory agree on the leading level but 
again have different corrections |69l |. 



V. THE ROLE OF DIELECTRIC JUMP AT CHARGED SURFACES 



So far we have assumed that the dielectric constant is uniform in space and equal for both solvent medium (where 
counterions are present) and charged surfaces. However, charged surfaces (macroions) usually have a dielectric con- 
stant, which is different from that of an aqueous solvent; for bio-soft materials, the dielectric constant, e', is usually 
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FIG. 7: a) Snapshots from MC simulations of two like-charged walls and counterions with (5^ — 0.95) and without = 0) 
a dielectric jump at the walls and for two different coupling parameters as indicated on the graph. The distance between the 
walls is A//i — 1.0. The figure shows the main simulation box, which is replicated periodically in lateral directions in the 
simulation model. The lateral box size is determined from the global electroneutrality condition as L/fi = \J nE,N, where A'^ is 
the number of counterions (here A'^ = 50). For the sake of representation, the extension of the system in lateral directions is 
rescaled in each case with L/^. 

b) The rescaled effective pressure between two like-charged walls as a function of the rescaled wall separation, A//i, as obtained 
from simulations in the presence of a dielectric jump of = 0.95 at both walls. Symbols correspond to different coupling 
parameters H — 0.5 (filled squares), 10 (open circles), 30 (filled triangles), 100 (filled circles), and connecting lines are guides 
to the eye. 

c) Global behavior of pressure between the walls shown in terms of the rescaled wall distance, A/fi, and the coupling parameter, 
H. Symbols are the zero pressure points; filled symbols show the stable bound state of the walls and open symbols show the 
meta-stable and unstable states 88]. Circles correspond to the case with a dielectric jump of 5^ = 0.95 at the walls, and squares 
represent the results for 6^ = (corresponding to Figure I^Js). 



smaller than that of water e « 80 {e.g. e' « 2 for hydrocarbon). This introduces a dielectric jump at charged 
boundaries, which can be treated theoretically using the method of image-charges. 

For a counterion of charge valency g at a charged wall of dielectric constant e', the image-charge is given by q' — qS^, 
where 

S. = (25) 

uniquely represents the dielectric jump. Intuitively, one expects that a dielectric jump of > leads to the depletion 
of counterions from the vicinity of the charged wall, since the image-charges have the same sign as counterions pll 
0, EH 13 ■ The depletion of counterions at a charged wall appears to be weak for small coupling parameters (the 
mean- field regime) and becomes significant for increasing coupling parameter (see below) [s^, 94]. The counterionic 
density profile for > shows a similar cross-over behavior as in the case of Se — (Section IIII C|l when coupling 
parameter is finite, i.e. it agrees with the strong-coupling prediction at small distances from the wall and follows 
the mean-field prediction for large distances (note that the strong-coupling theory in this case explicitly includes the 
image-charges). We shall not consider the case of a single charged wall here and only focus on Monte-Carlo results 
for the effective interaction between two like-charged walls in the presence of a dielectric jump [s^- 

Let us consider two like-charged walls of surface charge density —as and point-like neutralizing counterions confined 
in the space between the walls, which are located at a distance of A. We assume the same dielectric jump of (5^ > at 
both walls. In this case, one has to account for an infinite number of image-charges for each counterion, which leads 
to quite involved numerical calculations. We shall therefore proceed with an approximate description by taking into 
account only the first-order images (that is one image for each counterion in each wall). We show typical snapshots 
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from Monte-Carlo simulations of the two-wall system with and without a dielectric jump in Figure[7^. The depletion 
of counterions from the walls is clearly seen in these snapshots and appears to be stronger for larger coupling strength, 
since the repulsive interaction between counterions and their images grows with S. The effective pressure between 
the walls is shown in Figure [7|3 for a dielectric jump of ~ 0.95 (corresponding to water-hydrocarbon interface) 
for various coupling parameters. The image-charge interactions lead to a higher pressure at small wall separations 
as compared to the case with (5^ = (Figure EJi) . Also the range of distances at which the two walls attract each 
other is pushed to larger wall separations, which is more noticeable for larger coupling strength and implies a larger 
attraction at intermediate distances (see also Ref. 21J). However for sufficiently large wall separations, the dielectric 
effects weaken and the pressure becomes repulsive. 

The global behavior of the effective pressure is shown in Figure [7|; for = 0.95 (circles) and = (squares), 
where symbols show the points of zero pressure with filled symbols representing the stable bound state of the walls 
and open symbols representing the meta-stable or unstable states (corresponding to the local minimum and maximum 
of the free energy). As seen, the onset of attraction is shifted to somewhat larger coupling parameters (S « 30) in 
the presence of a dielectric jump, and also the bound state separation between the walls is larger and increases with 
the coupling parameter. However, one should be careful in drawing conclusions for the phase behavior of this system 
from the present data, since we have only considered the first-order images. The full numerical analysis of this system 
and comparison with extended strong-coupling theory (incorporating the dielectric jump) will be presented elsewhere 



VI. THE ROLE OF CURVATURE: CYLINDRICAL AND SPHERICAL MACROIONS 



In the preceding Sections, we considered planar systems, whereas in realistic situations, charged surfaces often have 
an intrinsic curvature. For simplicity, let us consider here only charged spherical and cylindrical macroions which are 
characterized by a single radius of curvature R. The radius of curvature sets a new length scale, which can introduce 
new features in terms of counterionic properties of the system. 

Intuitively, one may expect that when the radius of curvature of macroions is larger than the Gouy-Chapman 
length, R^ fi, the properties of the system remain qualitatively close to those of planar charge walls. The qualitative 
deviations from planar case may thus be expected for small R/ fi. As we shall see later, the geometrical symmetries 
of macroions, e.g. whether they be cylindrical or spherical, also play a role as they enforce the boundary conditions 
to which counterions are electrostatically coupled. Still it is useful to consider the dimensionless ratio of curvature 
radius, R, and the Gouy-Chapman length, /i, as 

which can characterize some important aspects associated with the curvature. For curved surfaces of uniform surface 
charge density CTs, we adopt the same definition for the Gouy-Chapman length as in Eq. Q, i.e. 

^ ( R/ (qiBT) charged cyhnders, 

M = t: — 7 = < (27) 

2TrqlB<Js y 2R^/{qtBZ) charged spheres, 

where r ~ 2'KasR is the linear charge density (in units of the elementary charge e) in the case of charged cylinders, 
and Z = AirasR^ is the total charge valency for charged spheres. It is important to note that in these cases, does 
not necessarily reflect the mean distance of counterions from the surface in contrast to the charged walls. 

For charged cylinders, the parameter ^ is quite well-known and is referred to as Manning parameter |89| | , which may 
be written as 

e = - = qiBT, (28) 

where we have used Eqs. (|26|l and (|27|l . By analogy we shall refer to the same ratio for charged spheres as Manning 
parameter, which reads 

Note that in a system with counterions of finite diameter, ad, the effective Gouy-Chapman length is larger than 
what one obtains for the same system with point-like counterions. Specifically, when counterions have a hard-core 
volume interaction with macroions, one has to use the hard-core radius of macroions in Eq. (|27l) . that is 

Rhc^R + fTci/2, (30) 
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which is larger than the actual radius, R, resulting in a reduced surface charge density for a given Z (or r) in the case 
of spheres (or cylinders). This leads to a larger Gouy-Chapman length, Eq. (|27|l . and a smaller coupling parameter, 
Eq. But the Manning parameter as defined in Eq. H26|l remains unchanged ^71] . 

A. Binding-unbinding transition of counterions 

The peculiar features emerging in the presence of charged curved surfaces are related to the behavior of countcrion 
at large distances from the surface. At equilibrium, counterions tend to diffuse away from macroions in order to 
maximize the entropy of the system, while at the same time, they are attracted energetically toward the macroion 
surfaces. 

For counterions at a charged sphere, the gain in entropy grows with the distance of counterions from the macroion 
center, r, qualitatively as ~ Inr for large distances. The energetic attraction, on the other hand, behaves like 1/r, 
which therefore is always weaker than the entropic repulsion experienced by counterions. Hence, counterions at a 
charged sphere tend to unbind completely and diffuse to infinity in the absence of confining boundaries. Thus the 
role of confinement becomes important in keeping counterions in the proximity of charged spheres. The confinement 
volume per sphere is related inversely to the concentration of spherical macroions in a solution. The qualitative 
considerations given above indicate that in the infinite- dilution limit (where the concentration of spheres tends to 
zero), the equilibrium counterionic density profile vanishes due to the complete de- condensation of counterions. Note 
that for a charged wall, the counterion-wall attraction grows linearly with distance, ~ z, overcoming the entropic 
contribution. The charged wall thus binds all its counterions, which is reflected by the fact that the density profile of 
countcrions-though extended to infinity as in the mean-field regime-is normalizablc to the total number of counterions 
(see Sections IIII Al and IIII B|l . We shall investigate these aspects further in the context of interaction between two 
spheres in Section IVlII 

The case of charged cylinders lies between the two cases of charged walls and spheres in that the energetic attraction 
of counterions to the cylinder grows logarithmically with the distance from the cylinder axis, r, i.e. in the same way 
as the entropic gain increases, ~ Inr. The competition between these two effects can result in a threshold binding- 
unbinding process in this geometry when the infinite-dilution limit is reached. In order to determine the threshold, 
one needs to consider the prefactor of both logarithmic contributions. The energetic attraction of a counterion to an 
infinitely long cylinder is given by /3U — 2q£BT\nr — 2^1nr (per fceT' = The entropic gain at large separations 

may be written as S/ftB — 21nr, since the cylindrical boundary implies a two-dimensional geometry. Comparing the 
two contributions, a threshold value of = 1 is obtained: for Manning parameter ^ > 1, the attraction wins and 
can lead to partial binding of counterions (with a finite density profile at the cylinder), whereas for ^ < 1, complete 
de-condensation of counterions is expected. This qualitative picture is actually supported by existing analytical 
results of mean- field [s^lsOl and strong-coupling 0| theory; both limiting theories give the same threshold of = 1. 
Recent numerical simulations [oJl show that this threshold is in fact universal and holds in all ranges of the coupling 
parameter S. Moreover, the threshold counterion-condensation process at charged cylinders exhibits a set of scaling 
relations, which are characterized by universal exponents j91|. 

In the case of two charged cylinders, the threshold process is expected to occur at Manning parameter ^ = 1/2 
as we shall discuss in Section rVIIII 

B. Attraction criteria for interacting spheres and cylinders 

The binding-unbinding behavior of counterions can drastically affect the effective interaction between macroions in 
solution, particularly, at low concentrations of macroions. When counterionic clouds at macroions become diluted due 
to the de-condensation process, counterion-mediated interactions are weakened and the effective interaction between 
macroions is dominated by their bare Coulombic repulsion. As mentioned above, the de-condensation process may 
in principle occur in all ranges of the coupling parameter, S, since it is regulated by the Manning parameter, ^, 
which is independent from the coupling parameter. Thus for curved surfaces-in contrast to planar systems-a large 
coupling parameter (S ^ 1) by itself does not necessarily indicate the regime of large electrostatic correlations, where 
strong-coupling attraction is expected between like charges. 

In order to specify the attraction-dominated regime for interacting spheres and cylinders, one can still employ 
a criterion similar to the Rouzina-Bloomfield criterion introduced for two charged walls (see below). But such a 
criterion should be supplemented by an additional condition on Manning parameter guaranteeing that a sufficiently 
large fraction of counterions condense in the vicinity of macroions. It is however difficult to establish this latter 
condition even for the simplest interesting cases of two spheres and two cylinders. Because it requires a detailed 
analysis of the binding-unbinding process of counterions in these systems, which is available only in the asymptotic 
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cases of mean field (S 0) [t^ 's^ and strong coupling (S oo) The mean-field theory is irrelevant for our 

purpose (as it does not include correlations), but as we shall see in Sections I V 111 and I V 1111 the strong-coupling theory 
can be used to obtain a quantitative prediction for the range of Manning parameters, where attraction may occur 
between macroions. In the following, we briefly mention the Rouzina-Bloomfield attraction criteria for spheres and 
cylinders assuming that the condition on Manning parameter is fulfilled. 



1. Like-charged spheres 

When counterions are highly condensed at spheres, we expect that correlation effects become dominant when the 
typical distance between counterions, a±, becomes much larger than the Gouy-Chapman length. Using the local 
electroneutrality condition, one may estimate a±_ as 

(31) 

(up to some numerical prefactor of the order unity), which also gives a measure of the correlation hole size around 
counterions at surface (FigurelS)). Note that in units of the Gouy-Chapman length, fi, we have aj_ = a^/ ji ^ (2S)^/^, 
where the coupling parameter for charged spheres may be written as 

" - g (-32) 

using Eqs. Q and (|27|l . The correlation- induced attraction is expected when the surface-to-surface distance between 
spheres, A, becomes smaller than the counterionic separation at the opposing surfaces, i.e. 

A < aj.. (33) 

In units of the Gouy-Chapman length, the above criterion may be written in terms of the coupling parameter as 
(A//i)2 < S, which is qualitatively similar to the attraction condition obtained for two char ged walls in Section 
in^^l The above attraction criterion was explicitly verified in simulations by AUahyarov et al. |30j and Linse et al. 
|32l l33| on like-charged spheres (Section [VII A|) . 



2. Like-charged cylinders 

For charged (parallel) cylinders, the attraction is similarly expected to arise for 

A < a^, (34) 

where A is the surface-to-surface distance of cylinders and is the typical separation between condensed counterions 
|4l| . In the case of highly-charged cylinders, counterions tend to accumulate in the intervening region between 
cylinders, where they line up on opposing surfaces along the cylinder axis forming a correlated inter-locking pattern 
pOj. The typical separation between counterions in this situation, a^, may be estimated from the local electroneutrality 
condition q = tQz, giving 

az = -, (35) 

T 

which agrees with the results obtained in recent simulations We shall discuss the application of criterion H34f) in 
Section IVIII Bl using numerical simulations of the two-cylinder system. 



VII. ATTRACTION BETWEEN LIKE-CHARGED SPHERES 



A. Numerical simulations 



Recently, there have been several simulations [2^129118(1 13^l33ll33 . l85j | investigating effective electrostatic attraction 
between like-charged spheres in the large coupling regime using multivalent counterions or in some cases, using low 
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Linse et al. \3^33j 
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Hribar et al. [34] 
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0.94 
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68 
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TABLE II: Parameters from simulations on highly-charged spheres: q is the charge valency of counterions with diameter a^i, 
Z is the charge valency of spheres with radius R, and ^b, M, C s-nd H are the Bjerrum length, Eq. Gouy-Chapman length, 
Eq. I27II . the Manning parameter, Eq. (1291 . and the coupling parameter, Eq. 1321 . respectively. L is the confinement box 
size and ^'^^ is the estimated attraction threshold discussed in Section IVll CI The last two columns show the equilibrium 
surface-to-surface distance obtained in these simulations, Asim (if explicitly measured), and the corresponding result from the 
strong-coupling theory. A*. Some of the numbers are given up to the order of magnitude, and the extracted values of A from 
simulations have a typical resolution of about lA. Note also that in estimating the values of E, fj , and a±, we account for the 
finite size of counterions by assuming that they have a hard-core interaction with macroions |2gl 13(1 13^ 133 . 13^ -see the note 
before Section IVra and Ref. 



dielectric constants |30| or considering the system at low temperatures • Like-charge attraction is reported in all 
these simulations for moderate to large coupling parameters (see Table . 

The strength of attractive force obtained between spheres is sufficie ntly lar ge t hat it can lead to closely-packed 
bound states (including large aggregates) between like-charged spheres |23, 123, Hfl HI, IH, yj, |33 ■ The bound-state 
corresponds to an attractive minimum in the potential of mean-force between spheres at small surface-to-surface 
separations 22i, 32^ Ji] • An interesting feature is that the attraction regime at small separations is separated by 
a pronounced potential barrier from a repulsion regime at large separations (see Fig. 1 in Ref. |28j|). On the other 
hand, the attractive minimum and the potential barrier are not robust and exhibit a dependence upon the size of the 
confinement volume 28] : for increasing confinement volume, the depth of the attractive minimum and at the same 
time the height of the potential barrier, decreases leading to a long-ranged repulsion between spheres in a sufficiently 
large confinement as expected. Though such a dependence on confinement volume appears to be quite weak (Section 

IviTTl . 

The existence of a potential barrier in the interaction potential of confined spheres can result in meta-stable bound 
states between two highly-coupled spheres 35] and also indicates a first-order phase transition (phase separation 
between a dilute and an aggregated phase) in the thermodynamic limit [2^ Is^, ■ 

We shall mainly focus on the attraction regime at small separations, where the attractive force between spheres 
maintains a typical equilibrium surface-to-surface distance of the order of the counterion diameter 2^ „29.. ,32., „33j . 
This equilibrium (or bound-state) separation may be obtained from the location of the minimum of the potential of 
mean force or equivalently from the pair distribution function of spheres. 

In Tabled we show the estimated values for the equilibrium surface-to-surface distance between spheres in these 
simulations, Agim, along with other parameters such as the coupling parameter, S, the Gouy-Chapman length, /x, 
as well as the estimated typical separation between counterions, a±, from Eq. H31|) . As seen the Gouy-Chapman 
length is quite small in these systems (/i ~ lA) compared with the counterion separation at spheres, aj_, resulting 
in a large coupling parameter. The strong-coupling attraction criterion H33I) is also fulfilled, since the equilibrium 
separation between spheres, Agim, is found to be smaller than a±. These observations suggest that the asymptotic 
strong-coupling theory is indeed relevant in the regime of parameters considered in these simulations. In the following, 
we shall present the strong-coupling predictions and compare them with the simulation results. 



B. Effective sphere-sphere interaction: Asymptotic strong-coupling theory 

Let us consider a system of two like-charged spheres (Figure ^ that are in general confined in a cubic box of 
edge size L inside which the global electroneutrality condition is satisfied, i.e. Nq = 2Z, where N is the number of 
counterions and Z = AirR^as is the (absolute value of the) charge valency of each sphere. We assume that counterions 
(of diameter ad) have also a hard-core excluded-volume interaction with the spheres. 

We are interested in the limit of large couplings. As mentioned in Section llV Bl the leading contribution to the free 
energy for S ^ oo involves only the one-particle contributions. The strong-coupling free energy for the two-sphere 
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FIG. 8: Two identical like-charged spheres of radius R axe considered at center-to-center distance of D in a cubic box of edge 
size L. The bare charge of spheres with uniform surface charge density —as (in units of the elementary charge, e, and assumed 
to be negative) is compensated by the total charge of counterions of charge valency g (and diameter ad). 



system may be written as (up to an irrelevant additive term) [tH 

P^SC , /" 1 , , -pu{x,y,z) 
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In / dx dy dz e~ 
'v 



(36) 



where D = is the center-to-center distance of spheres in units of the Gouy-Chapman length. The first term 

in Eq. (|36|l is the bare Coulombic repulsion of spheres with Manning parameter ^. In the second term, /3m is the 
single-particle interaction energy of counterions with the spheres. 
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in which fi = ri/fj, and r2 — ^^e the rescaled radial distances from the centers of the two spheres labeled 1 

and 2 (we choose the frame of reference in the middle of the box such that ri^2 = [{x ± I?/2)^ + + z^]-^/^). The 
spatial integral in Eq. H3()|l runs over the volume accessible for counterions, i.e. inside the cubic box excluding the 
two spheres (note that a shell of thickness crci/2 around each sphere is also excluded due to the hard-core interaction 
between spheres and counterions). 

The energetic and entropic contributions from counterions enter on the leading order through the second term in Eq. 
(|36|l . It is important to note that the main qualitative features regarding the binding- unbinding behavior of counterions 
is reproduced by this term. In particular, for very large confining box L — > cx), the single-particle partition function 
Zi(Z?,^, L) — Jy d'^r exp(—Pu) diverges with the box volume as Zi V = for any given Manning parameter 
(since the integrand is always positive and bigger than one). Thus the distribution of counterions around the spheres, 
~ exp(—f3u)/Zi, as well as the component of the force contributed by counterions, ~ dlnZi/dD ^ L~^, vanish in 
the limit L — > 00. This leads to a pure repulsion between unconfined spheres as expected. Note that this repulsion 
regime does not represent the strong-coupling situation associated with large electrostatic correlations. The fact 
that both the de-condensation process and the repulsive regime for unconfined spheres are consistently captured by 
the asymptotic strong-coupling free energy indicates that the 1/S-expansion scheme used to obtain the asymptotic 
contribution for S — )■ cxd 69] is not only based on energetic considerations, and can also account for entropic effects 
on the leading order [691, (tJ . 

Now let us consider the spheres in a finite confinement volume and investigate the strong-coupling prediction 
for the interaction free energy (|36|l . For small Manning parameters, the free energy exhibits only a long-range 
repulsion (Figure |5^), but as Manning parameter exceeds a threshold value of a local minimum is developed 
at small separations indicating a short-range attraction and a meta-stable bound state. As seen (Figures I^Jd and c), 
this attraction regime is separated from the large-distance repulsion regime by a pronounced potential barrier. For 
increasing Manning parameter, the attractive local minimum becomes deeper than the large-distance minimum, and 

the potential barrier disappears beyond a second threshold of ^c^''-see Figure The preceding features indicate 
a discontinuous unbinding transition between a closely-packed bound-state and a repulsion-dominated state of two 
like-charged spheres by varying the Manning parameter. The generic form of the free energy also agrees qualitatively 
with numerical findings |22j as discussed in Section IVlI Al 
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FIG. 9: The strong-coupling free energy of the two-sphere system, Eq. (1361 . plotted as a function of the rescaled center-to- 
center distance, D / ^, for Manning parameters a) ^ = 1.5, b) ^ = 3.4, c) ^ = 3.6 and d) ^ = 4.0 (the box size is L/ ji — 100). 
The location of the local minimum at small separations is marked by an arrow. A potential barrier is found in the range of 
Manning parameter ^i^\L) < £, < where we have ^c^' « 3.3 and ^c^' ~ 3.8 for L/ii — 100. The values of the two free 

energy minima become equal at ^ ~ 3.5. 




A digram representing different regimes of attraction and repulsion is shown in Figure 1101 where we have plotted 
the locations of the minima (solid curves) of the strong-coupling free energy and also the location of its maximum 
(dashed curves) as a function of Manning parameter, ^. (The locations of the two threshold Manning parameters, 

S}c^^ and are shown by arrows.) One can show that for ^ <C 1, the location of the repulsion-dominated minimum 
at large separations (Figure |5^), D*, scales linearly with the box size as 

(38) 

when the box size tends to infinity L ^ oo |23- On the other hand, for large Manning parameter ^ ^ 1, the location 
of the attraction-dominated minimum at small separations (Figure Eli), saturates to a value independent from 
the box size; in this case, the equilibrium surface-to-surface distance of spheres. A* = I?* — 2i?, follows from Eq. 
approximately as 

A, eeD, -2i?waci + yA' + C'(M^), (39) 

where is the counterion diameter. This minimum corresponds to a highly-condensed state of counterions in the 
intervening region between spheres jTH . The attractive force induced between spheres in this regime (^ 3> 1) may also 
be calculated from the strong-coupling free energy (|36() : at small separations {D ^ 2i?), the force takes the following 
analytical form 

This limiting attractive force (for S ^ 1 and ^ 3> 1) is independent from the temperature and results only from 
energetic contributions |71j. The expression l|4()|l qualitatively agrees with the results obtained by Shklovskii using 
the Wigner-crystal model |65l |. 

The predictions of the strong-coupling theory for the bound-state separation of two attracting spheres. A*, has been 
compared with numerical simulations in Table m (note that here A^* has been calculated by numerical evaluation of 
the free energy H36(l for the corresponding simulation parameters. The analytical expression 1)39(1 may also be used, but 
it gives an approximate value up to the first order in jx). As seen there is a reasonable semi-quantitative agreement 
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FIG. 10: Regimes of attraction and repulsion in the system of two like-charged spheres and counterions as obtained from the 
asymptotic strong-couphng theory. The solid curves show the rescaled equilibrium center-to-center distance between spheres 
as a function of the single-sphere Manning parameter, ^ (Eq. 12911 ). for a confining box of rescaled size L/^i = 100. The 
dashed curve corresponds to the maximum of the strong-coupling free energy and the dotted vertical line shows the Manning 
parameter for which the values of the two minima of the free energy are equal (see Figure |3J . The locations of threshold 
Manning parameters ^'^^ and ^c^' are shown by arrows. The dot-dashed line shows the contact separation D = 2R. Inset: 
The attraction threshold for two like-charged spheres, £}c \ increases weakly with the rescaled box size, L/^i. It exhibits a 
logarithmic dependence according to Eq. II42II . which is shown by the solid curve. 



between the theoretical predictions and the simulation results. As mentioned in Section IVll Al the equilibrium 
surface-to-surface distance in these simulations appears to be about the counterion diameter (see the Discussion in 
Ref. 113 )• This also follows from the strong-coupling prediction, Eq. since for highly-charged spheres, the 

Gouy-Chapman length is in fact small compared with the counterion diameter. One should also note that the lateral 
separation of counterions at spheres in the simulations is typically larger than the counterion diameter indicating 
that the excluded- volume interaction between counterions is not a dominant effect. The volume interactions between 
counterions enter only in the hi gher - order corrections to the asymptotic theory and can lead to additional attractive 
components between macroions |30l |40| . 



C. Attraction threshold 

The strong-coupling theory also allows to obtain an analytical estimate for the regime of Manning parameters where 
attraction is expected between spheres. This regime may be specified by 

c>e^ (41) 

where the threshold Manning parameter ^c^"* actually depends on the confinement size, s}c'^ — s}c'\l)- As shown in 
the inset of Figure s}c^^ [L) increases almost logarithmically with the box size as 

fW(L)«a + Mn(^^), (42) 

where a ~ 0.55 and b ~ 0.6 are obtained by fitting to numerically-determined SC predictions (symbols). The estimated 
values of ^c^^ are shown in Table ^ for the given simulation parameters, which show that these simulations indeed 
exhibit the attraction regime (|41|l . The weak dependence of the attraction threshold, on the confinement size can 
also explain the stability of compact clusters of spheres in quite large confinements (i//i ^ 1) 28, 2^ l3(ill8^ l33. 34, 3^, 
since the Manning parameter only needs to exceed a moderate value (^ In L/ /i) for like-charged spheres to fall into 
the attraction-dominated regime. 
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FIG. 11: Two identical and parallel charged cylinders of radius R are considered at axial separation of D and in a square box 
of lateral edge size L; cylinders have a length of H, which is assumed to be infinitely large. The bare charge of cylinders with 
uniform surface charge density —as (in units of the elementary charge, e, and assumed to be negative) is compensated by the 
total charge of counterions of charge valency q (and diameter CTci). 



VIII. ATTRACTION BETWEEN LIKE-CHARGED CYLINDERS 



A. Effective cylinder-cylinder interaction: Asymptotic strong-coupling theory 

In this Section, we shall consider the interaction between two like-charged cylinders in the limit of large coupling 
parameter S ^ oo. For simplicity, we assume that the cylinders are infinitely long (with length H), and that they 
are confined in a square box of edge size L (see Figure [TT|) . The electroneutrality condition holds inside the confining 
box, thus we have qN = 2tH, where N is the number of counterions and r — 2'KasR is the linear charge density of 
each cylinder. We assume that counterions and cylinders have also a hard-core excluded-volume interaction. 

The strong-coupling free energy of this system (up to an irrelevant additive term) follows as 0| 

f^^^ ^ -^\nD -\n J dxdye-^"(^'2^), (43) 

where D ~ D / j^i is the axial separation of the cylinders in units of the Gouy-Chapman length. The first term in 
Eq. H43|l is the bare repulsion of cylinders and in the second term, fiu is the single-particle interaction energy of 
counterions with the cylinders, 

/3u = 2^(lnfi +lnf2), (44) 

where fi — r\l and f2 = ^2! [i- are the radial distances from the axes of the two cylinders (in the a;?;-plane perpendicu- 
lar to the cylinders axes)-we choose the frame of reference in the middle of the box such that r\^2 — [(a;±£'/2)^-)-?/^]^/^. 
The integral in Eq. H43|l runs over the volume accessible for counterions inside the box excluding the two cylinders 
(and a shell of thickness crci/2 around each cylinder corresponding to the closest approach distance of counterions). 

As in the case of two spheres (Section IVIIBII . the effective interaction between like-charged cylinders is influenced 
by the binding-unbinding behavior of counterions incorporated on the leading order in the single-particle partition 
function Z\ {D, ^, L) = Jy d^r exp(— /3m) in the second term of the free energy (|43|l . For further analysis of this behavior 
in the two-cylinder system, let us first consider the limit of very large box size L 00. In this limit, Zi scales with 
the box size as Zi ~ which may be seen simply by rescaling the spatial coordinates with i as a; — > x/L, etc. 

Thus for ^ < 1/2, Zi diverges and consequently, the distribution function of counterions, exp{—(3u)/Zi, vanishes 
indicating de-condensation of counterions from the two cylinders. The counterion-mediated force between cylinders, 
~ dlnZi/dD, tends to zero as well; thus the cylinders only repel each other in the limit L 00. In contrast, for 
Manning parameter ^ > 1/2, counterions bind to the two unconfined cylinders with a finite density profile, and thus 
can produce attraction between them. The condensation threshold obtained above for two cylinders, = 1/2, agrees 
with the classical result due to Manning [t^ quantitatively. 

An interesting question is to obtain the regime of Manning parameters, where the attraction emerges between 
cylinders. Intuitively, we expect that the effective attraction sets in somewhat above the counterion-condensation 
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Dl\x. Dl\x. 

FIG. 12: The strong-coupling free energy of the two-cyhnder system as a function of the rescaled axial distance, -D//i, for 
Manning parameters a) ^ = 1.0 and b) ^ = 0.3 and for several box sizes as indicated on the graph. Arrows show the 
approximate location of the minimum of the free energy, which for large Manning parameters is nearly independent of the box 
size and reflects a closely-packed bound state. For small Manning parameter ^ < = 2/3, the location of the minimum tends 
to infinity with the box size reflecting a continuous unbinding transition for like-charged cylinders. 



threshold = 1/2. Because right at the condensation threshold, there is an unbalanced bare repulsion between the 
cylinders, and a finite fraction of condensed counterions is needed to compensate this repulsion. It follows from a 
more detailed analysis of the strong-coupling free energy (|43f) that the attraction actually emerges when [tJ 

C> = 2/3, (45) 

which gives a universal threshold for attraction between unconfined cylinders. 

The attraction threshold between two like-charged cylinders has also been investigated using other methods. Anal- 
ysis of Ray and Manning |78j based on the classical counterion-condensation theory [s^ predicts attraction between 
two like-charged cylinders for ^ > 1/2, which coincides with the onset of counterion condensation. It should be noted, 
however, that the attraction mechanism involved in their theory is not based on electrostatic correlations, but features 
a mean-field covalence-like binding process. Arenzon et aL's study 113,123 based on a structural-correlations theory 
(which also accounts for counterion condensation) predicts attraction for ^ > 2. Numerical simulations |27l liol IllL 14^ . 
on the other hand, give attraction for the range of Manning parameters ^ > 0.8, but have not yet specified the at- 
traction threshold precisely fSection lVIII Bp . 

For two cylinders in a finite confinement volume, one can determine the attraction and repulsion regimes by 
evaluating the strong-coupling free energy Eq. H43|l . The typical form of the free energy H43|l is shown in Figure 
El for both large (C > ■Cc = 2/3) and small (^ < = 2/3) Manning parameters. For large Manning parameter 
(Figure EJl), the free energy exhibits a long-range attraction and a local minimum at small separations, which is 
quite insensitive to the confinement size as counterions are localized mostly in the proximity of the cylinders in this 
regime. The analytical form of the attraction force (per unit length of the cylinders H) may be estimated for large 
Manning parameters ^ 3> 1 [tH : in the zero-temperature limit (^ — *■ oo), we have: 



F{D) 
H 




(46) 

D w 2i?, 



for large and small axial separations respectively. This limiting attractive force (which is obtained for 5 1 and 
^ ^ 1) is independent from the temperature and originates from a purely electrostatic origin. 

For small Manning parameter ^ < = 2/3 (Figure EId), the interaction free energy of cylinders is dominated 
by their bare repulsion as counterionic clouds become increasingly diluted around the cylinders. The location of the 
minimum of the free energy for ^ <C 1, -D*, tends to infinity with the box size as « L/^/tt. 

In brief, one may specify the attraction and repulsion regimes of two like-charged cylinders by considering the 
location of the minimum of the free energy (|43|) as shown in Figure 113b for several different box sizes (the region 
below each curve shows the repulsion regime and above that is the attraction regime). For ^ S> 1, the equilibrium 
surface-to-surface separation of the cylinders. A*, may be obtained approximately as 

A* =D*-2i?«aci + ^M + C'(M^), (47) 
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FIG. 13: a) Regimes of attraction and repulsion in the system of two like-charged cylinders and counterions as obtained from 
the asymptotic strong-coupling theory. The curves show the rescaled equilibrium axial separation between cylinders, -D*//i, as 
a function the single-cylinder Manning parameter, ^ (Eq. for confining boxes of different size as indicated on the graph. 

The curves are obtained numerically by minimization of the free energy 14311 . In an unconfined box, the cylinders continuously 
unbind at the threshold Manning parameter = 2/3 (shown by an arrow). The thick gray curve corresponds to D = 2R, 
where the two cylinders are at contact. 

b) The rescaled equilibrium surface-to-surface separation of the cylinders, A,/fi = (D« — 2R)/fi, as a function of Manning 
parameter, 5. Symbols are simulation data for the Rouzina-Bloomfield parameter (from top): 7rb = 3 (crosses), 10 (triangle- 
downs), 15 (triangle-ups) , 30 (stars), 40 (squares), 50 (open diamonds) and 60 (filled diamonds) ISjl . The solid curve is the 
strong-coupling prediction obtained from Eq. 11431 . Dashed curves are guides to the eye. 



where ad is the counterion diameter. For typical strongly-coupled systems (with small Gouy-Chapman length), the 
equilibrium surface separation is predicted from Eq. H47fl to be of the order of the counterion diameter. 

When Manning parameter decreases down to the attraction threshold = 2/3, the two cylinders unbind from 
each other for L — s- cx) (solid curve in Figure IT^ ). But in contrast to like-charged spheres, the unbinding transition 
of cylinders occurs in a continuous fashion exhibiting a universal scaling exponent for the diverging axial distance of 
cylinders ^ 

D^^i^-^c)-", (48) 

where the exponent is found as 

a = 3/2. (49) 



B. Comparison with numerical simulations 

Several numerical simulations have been reported on char ged cylinders [27L l40l l4ll |4^ and also on more detailed 
models, which incorporate the charge pattern of DNA [26l l36l|. Attraction is found in a wide range of Manning 
parameters (including ^ « 1.0) and for moderate to large coupling parameters. 

Gr0nbech- Jensen et al. 27] showed that the effective interaction between two parallel like-charged cylinders in 
the presence of divalent counterions exhibits repulsion at very small distances close to contact and attraction at 
intermediate distances by means of which a closely-packed bound state is maintained at small surface separations. At 
large distances the effective interaction turns to repulsion. Similar results have been reported in simulations by Lee 
et al. (42.] . In these simulations, charges on cylinders and counterions interact with Coulombic interactions as well 
as short-range excluded- volume interactions. Recent simulations by Deserno et. al [i^ indicate an interplay between 
electrostatic and excluded-volume interactions, which may drastically influence the effective electrostatic interaction 
between like-charged cylinders. Specially for large linear charge densities (when charge separation on the cylinder is 
smaller than the counterion size), the attractive electrostatic force between cylinders weakens and eventually turns to 
repulsion for sufficiently large coupling strength ^3 • In this case, the total effective force between the cylinders is still 
attractive due to an attractive component from excluded-volume interactions. This behavior is not yet completely 
understood and appears to be different from a simple depletion mechanism 1^. A similar situation may be present 
in the system studied in Ref. [27i] . however, the detailed analysis of the force components is not reported. 
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Intuitively, we expect that excluded-volume effects become important when typical distance between counterions 
on cylinders surface, = q/T fSection lVllj 2(1 . is smaller than the counterion diameter, Gd, i-e. for < a^i- In this 
case, excess accumulation of counterions in the intervening re gion between cylinders, which is favored energetically 
and leads to the strong-coupling attraction 7)] , is prohibited |4Cl'| . In order to investigate the electrostatic features 
of like-charge attraction in the two-cylinder system and to avoid complications arising from volume interactions, we 
consider recent Molecular Dynamics simulations performed by Arnold and Holm, which exclude volume interactions 
between counterions though counterions still retain a soft-core excluded-volume repulsion with the cylinders |4lj |. 
This is important for the forthcoming comparison made with the strong-coupling predictions, since the leading-order 
results at large couplings only involve counterion-cylinder interactions f Section IVIII|I . The simulation model has a 
geometry similar to what we have shown in Figure 1111 where periodic boundary conditions in z direction is used. 
(For convenience the simulated model employs a cylindrical outer boundary of diameter 8D. For final comparison, 
the theoretical curves are also calculated using a similar constraint, though for simplicity, a square box of edge 
size L = 8D is used, which as explicitly checked does not influence the following results in the considered range 
of Manning parameters The soft-core repulsion in the simulations is also chosen strong enough such that it 

effectively prevents counterions from penetrating the cylinders, which makes the simulation model comparable to the 
theoretical model with a hard-core volume interaction. In the simulations, cylinders are kept at fixed surface-to-surface 
distance, A, with fixed linear charge density, r, and counterion valency, q, but the Bjerrum length and the cylinders 
radius, R, are varied. Hence, the Gouy-Chapman length, /i (Eq. (|27ll '). varies accordingly allowing to span the phase 
space of the system and specify the attraction and repulsion regimes. 

In Figure [CTj. we show the simulation results for the equilibrium surface-to-surface distance of cylinders in rescaled 
units, A,//i, as a function of Manning parameter ^ = qiBT (Eq- (EHJ)- Different symbols correspond to different 
values of the Rouzina-Bloomfield parameter defined as 

7RB = (50) 

which is fixed for each simulation data set. This dimensionless ratio gives a measure of the deviations from the 
strong-coupling regime, since it represents the ratio between the estimated correlation hole size, = q/r |40l|. and 
the surface-to-surface distance between cylinders (see Section FVI B 2\i 92]. Thus according to Rouzina-Bloomfield 
criterion, large 7rb = ciz/ A. corresponds to a strongly-coupled system with highly correlated structures on opposite 
surfaces of cylinders that may give rise to attraction between them. This behavior is clearly supported by the present 
simulations as seen in Figure ITSb . For increasing 7rb (from about 3 up to 60), the equilibrium separation between 
cylinders decreases and tends to the strong-coupling prediction indicating a closely-packed bound state at large 
Manning parameters ^ > 1, where the equilibrium surface-to-surface distance is of the order of the Gouy-Chapman 
length. 

The quantitative agreement with the strong-coupling prediction is obtained for the whole range of Manning pa- 
rameters studied in the simulations. Due to convergence limitations, the simulations so far have been limited to the 
range of Manning parameters ^ > 0.8. It becomes more difficult to obtain good data for smaller Manning parameters 
as the distance between cylinders rapidly increases. Nonetheless, the excellent convergence of the present data to the 
strong-couphng curve suggests an attraction threshold of about = 2/3 obtained in Section IV 111 Al 

IX. CONCLUDING REMARKS AND DISCUSSION 

We have discussed the effective electrostatic interaction between like-charged macroions in the regime of large 
coupling parameter, S, which is achieved for large counterion valency, large charge densities on macroions, at low 
dielectric constants or low temperatures. In this regime, interactions between macroions are dominated by strong 
counterionic correlations: counterions form highly correlated layers at macroions, including Wigner crystals for suf- 
ficiently large coupling, that lead to attractive forces between like-charged macroions at sufficiently small surface 
separations. The relevant length scale is set by the typical distance between counterions at macroion surfaces (the 
correlation hole size) below which like-charged surfaces couple to each other electrostatically. This energetic coupling 
is essentially mediated by single counterions sandwiched between opposite surfaces of macroions. At larger surface 
separations, however, correlation effects and the attraction strength are reduced and multi-counterion interactions 
play a significant role. At very large separation (as long as the coupling parameter is finite), the interaction between 
macroions is dominated by mean-field features and eventually turns to repulsion. 

The strong-coupling attraction appears with the same mechanism for planar and curved surfaces. For curved 
surfaces, however, attraction is found only when a sufficiently large number of counterions are condensed near the 
surface; the condensation process is regulated by the ratio between the radius of curvature, R, and the Gouy-Chapman 
length, /i, that is the Manning parameter ^ = i?//Lt. Therefore, the attraction regime for charged cylinders and spheres 
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is characterized by both a large coupling parameter S, which is specified via the Rouzina-Bloomfield criterion, and also 
by a large Manning parameter, which has to be determined by considering the condensation process of counterions 
in the given system of macroions (Sections IvTIvTTTIi . 

We reviewed recent analytical results obtained in the limit of S — > c», which gives rise to the asymptotic strong- 
coupling theory. This theory incorporates single-particle contributions on a systematic level and may be extended 
to finite coupling regimes by considering higher-order corrections via a 1/S-expansion series (virial expansion) |69| |. 
The predictions of the strong-coupling theory were considered for the effective interaction between two charged walls 
and also for two spherical and cylindrical macroions. This asymptotic theory yields a long-ranged attractive force 
between like-charged objects that is of constant strength for two walls, but varies with distance, D, as ~ 1/D for two 
cylinders and as ^ l/D"^ for two spheres. Note that the leading-order free energy contains in addition to energetic 
contributions also (repulsive) entropic contributions from counterions (see Eqs. (I15f) . (|36|l and (|43l) ). which vanish 
at zero temperature. The zero temperature limit is in fact contained within the strong-coupling theory and may 
be obtained by taking the limit /x — > (or ^ — > oo for curved surfaces), where the free energy reduces to a purely 
energetic contribution giving a long-ranged attractive force as specified above (see Eqs. (fT5|l . (BHll . and Ref. [71]). 
As we know from analytical and numerical studies of planar walls 39, '69], the range of strong-coupling attraction 
is reduced substantially at finite couplings due to higher-order corrections. (The higher-order corrections have not 
yet been calculated for the system of two cylinders and two spheres.) By comparing the strong coupling predictions 
with recent simulations, it was shown that the attraction regime and the closely-packed bound state of macroions 
can be described by this limiting theory on a quantitative level |37, 38, 39, 41, 71j. Note that the agreement with 
simulations is limited to the regime of surface separations determined by the Rouzina-Bloomfield criterion. This 
regime particularly includes typical coupling parameters of S ~ 10^ and Manning parameters of ^ 10 (for spheres) 
and ^ > 1 (for cylinders), which are accessible in usual experimental systems (see Tables III and Hill . 

The asymptotic strong-coupling theory was also considered to estimate the threshold Manning parameter above 
which the attraction is expected for two like-charged spheres and cylinders. The attraction threshold is captured within 
the asymptotic theory because it accounts for the entropy-driven de-condensation process of counterions on the leading 
order. One should note, however, that the de-condensation regime does not involve electrostatic correlations and thus 
in general, the strong-coupling theory may be only qualitatively valid when counterions de-condense at low Manning 
parameters. Yet strong-coupling predictions were shown to remain in a reasonable agreement with simulations on 
two cylinders for decreasing Manning parameter close to de-condensation threshold [4l|. (Moreover, the predicted 
de-condensation threshold itself agrees quantitatively with the standard Manning results for one and two cylinders 
[zIjIs^-) Further numerical and analytical studies are useful to specify the validity of the strong-coupling predictions 
at low Manning parameters, including the predicted unbinding behavior of cylinders and spheres. 

We did not discuss possible thermodynamic phase transitions triggered by attractive forces in the system of charged 
plates, spheres and cylinders. (Note that the binding-unbinding behaviors discussed in Sections IVIII and lVIIII reDresent 
thermodynamic phase transitions only in the case of infinitely long cylinders; for spheres (or short cylinders) there 
will be additional entropic contributions to the free energy from the sphere-sphere distance coordinate, which are not 
considered within the present model.) A first-order unbinding transition was predicted to occur in the system of two 
like-charged walls at the coupling parameter S w 17 that has also been compared with experimental observations |39j |. 
Also there has been indication of an attraction-induced phase separation in the system of like-charged spheres from 
recent numerical simulations |27l Is^ Is^ . The systematic study of such phase transitions still remains a challenging 
subject (see Ref. and references therein). 

Another interesting problem is to examine the influence of additional salt on the interactions in the strong-coupling 
limit. The present results are expected to remain valid at sufficiently small salt concentrations (large screening 
length). Qualitatively, one can associate the size of the confinement box considered in the present models with 
the screening length. Thus, addition of salt is expected to matter particularly close to the unbinding threshold 
of macroions. Other relevant subjects include the effect of finite polymer stiffness, the discrete charge pattern of 
macroions [2^ [Sa, |53, IH, |63, IS^, |63, l^^ , and bundling of many charged polymers [13, [la 1^ IH> UM i'^ the strong 
coupling regime, which constitute interesting applications for the future. 
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